Weak Lensing of the CMB: Sampling Errors on B-modes 
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The B modes generated by the lensing of CMB polarization are a primary target for the upcoming 
generation of experiments and can potentially constrain quantities such as the neutrino mass and 
dark energy equation of state. The net sample variance on the small scale B modes out to I = 2000 
exceeds Gaussian expectations by a factor of 10 reflecting the variance of the larger scale lenses 
that generate them. It manifests itself as highly correlated band powers with correlation coefficients 
approaching 70% for wide bands of Al/l ~ 0.25. It will double the total variance for experiments 
that achieve a sensitivity of approximately 4 /iK-arcmin and a beam of several arcminutes or better 
This non-Gaussianity must be taken into account in the analysis of experiments that go beyond first 
detection. 



I. INTRODUCTION 

As a step on the road toward the ultimate goal of de- 
tecting primordial gravitational waves, upcoming cosmic 
microwave background (CMB) polarization experiments 
will target the distortion to the acoustic polarization in- 
duced by gravitational lensing. As with the polarization 
induced by gravitational waves, the gravitationally lensed 
polarization contains a component with handedness, the 
so-called B mode component 0. Unlike gravitational 
waves, gravitational lensing provides a guaranteed sig- 
nal. In the standard cosmological model, the predicted 
amplitude of the B modes can only vary at the tens of 
percent level within current constraints Moreover 
these fine variations provide an opportunity to measure 
the dark side of the universe, namely the dark energy 
and neutrino dependent growth of structure, as well as 
another handle on the reionization optical depth 

Although both the intrinsic distribution and the den- 
sity perturbations that lens the CMB are expected to be 
Gaussian, the lensed distribution is non-Gaussian at sec- 
ond order in the perturbations. The non-Gaussianity is 
therefore relatively small in the temperature distribution 
m. However because the B modes are generated by the 
lensing itself, its non-Gaussianity is a first order effect 
but fortunately one that is precisely calculable. Gravi- 
tational lensing therefore also provides a unique testing 
ground for experimentally extracting a non-Gaussian sig- 
nal in the presence of foregrounds and systematic errors. 
Ultimately, the non-Gaussianity of the lensed polariza- 
tion also provides the key to mapping the dark matter 
|6|, i3 and hence the separation of the lensing and gravi- 
tational wave B mode components 0, . 

For the upcoming generation of experiments, the non- 
Gaussianity will provide an important source of uncer- 
tainty for power spectrum measurements. Fisher infor- 
mation studies have shown that the information con- 
tained on cosmological parameters in the B mode power 
spectrum under the Gaussian approximation unphysi- 
cally exceeds that contained in the two underlying Gaus- 



sian fields 0. 

In this paper, we study the origin and quantify the 
impact of the non-Gaussian sample variance on B mode 
power spectra measurements. The basic reason for the 
large sample variance is that the fluctuations that lens 
the CMB are mainly on degree scales. All of the ar- 
cminute scale B modes fluctuate jointly with the lens 
and so precision measurements will require many degree 
scale patches not simply many arcminute scale patches. 

We begin in [^by briefly reviewing the generation of 
B modes through gravitational lensing. We calculate the 
non-Gaussian sample covariance in illlll and explore its 
impact on measurements in ^IVI We conclude in 
For illustrative purposes we employ throughout a fidu- 
cial cosmology that is consistent with WMAP determi- 
nations: an initial scale invariant spectrum of curvature 
fiuctuations with amplitude 6,^ = 5.07 x 10~^ {as = 0.91, 
T ~ 0.17), a baryon density Ui,h^ = 0.024 and a matter 
density ^mh^ = 0.14 in a flat J^a = 0.73 cosmology. 



II. B MODES 

Weak lensing by the large-scale structure of the uni- 
verse remaps the polarization field or equivalently the 
dimensionless Stokes parameters Q(n) and U{n) as P, 

mm 



[Q±iU]{n) = [g±z(7](n + V0(n)), 



(1) 



where n is the direction on the sky, tildes denote the un- 
lensed field, and the defiection angle Vt/i is the gradient of 
the line of sight projection of the gravitational potential 
*(x,Z?), 



(n) ^-2 dD 



{ Ds - D) 



^{Dn,D), 



(2) 



where D is the comoving distance along the line of sight in 
the assumed flat cosmology and Dg denotes the distance 
to the last-scattering surface. In the flat sky approxima- 
tion, the Stokes parameters can be decomposed into E 
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FIG. 1: B power spectrum = Z^Cf ^/27r and the diflFeren- 
tial contribution to the total B variance Atot from power in 
the lensing field (f> with wavenumber L [see Eqn. ©]. Even 
though the B power peaks at / ~ 10"^, it gains half its contri- 
bution from power in the lensing field at L < 450. Sampling 
errors in B therefore reflect the larger scale variations in <j). 



and B Fourier modes as 



[Q±iU]{n) = - 



(2^)2 

where 1 = (/ cos (/S/, Z sin (/?;) and likewise 



(n) = 



(27r) 



(/.(l)e^'-", 



(4) 



Even in the absence of an intrinsic B, lensing will gen- 
erate _B as 



B[\) 



2U 



;E(\'m~\')W{\,\'), (5) 



where we have expanded Eqn. in the leading order 
gradient approximation and the mode coupling weight is 



W(l,l') = [l'-(l-l')]sm2((^/-</7r) 



(6) 



Upcoming experiments will mainly measure the power 
spectrum of the B modes. In the flat sky approxima- 
tion, the power spectra of a statistically homogeneous 
field X{\) is given by 



{X*{\)X{\')) = (27r)2(5(l - l')Ci , 



(7) 



and statistical isotropy requires C\ = Ci. It follows that 
(e.g. [11) 
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FIG. 2; Fractional deviation of band powers from the ensem- 
ble average for three Monte Carlo realizations (lines) com- 
pared with the expected rms deviation (shaded region) for a 
Gaussian field of the same power spectrum. The realizations 
show correlated deviations in the recovered band powers at 
high I. 



given that the intrinsic polarization from z ^ 1000 is 
essentially uncorrelated with the lensing potential from 
z < 3. The B power spectrum is essentially a convolution 
of the E and <f) power spectra. 

The total power in the B field arises from power in the 
6 field as 



^,Cr = (0.46mK/T)^ 



(9) 



(2^)2 



where the value in parentheses is the rms in the fiducial 
cosmology and T = 2.725 x lO^^K is the CMB tempera- 
ture. For reference, the rms of the I < 2000 low-pass fil- 
tered B field is 0.43/iK. We will typically take this value 
as the maximum I estimated for illustration purposes. 

In Fig.^ we show the B power spectrum in the fiducial 
cosmology. Note that the power spectrum peaks at I ~ 
10'^ reflecting the power in the underlying unlensed E 
power spectrum. However this mode coupling in Eqn. (jS]) 
is achieved through power in the potential C^"^ across a 
broad range in L. The B field acquires half its power 
for L < 450. In other words, degree-scale gravitational 
lenses give rise to the B power at the 10' scale. 

This fact is the key to understanding the non-Gaussian 
covariance. In a given degree scale patch of sky there is ef- 
fectively only a single lens contributing to the small scale 
B power. The amplitude of all of these smaller scale B 
modes will covary with the amplitude of the larger scale 
lens. Hence reducing the sampling errors for arcminute 
scale B modes requires many independent degree scale 
patches of sky. If the B field were Gaussian, such a 
constraint would require only an equivalent number of 
arcminute scale patches of sky. 
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We illustrate this effect with Monte Carlo realizations 
of the polarization field and lenses from the power spec- 
tra of the fiducial cosmology. The Monte Carlo simu- 
lations were performed using a square patch of sky of 
side length 22.9 degrees, and with fields sampled on a 
grid with 1.3 arcminute spacing. Periodic boundary con- 
ditions were used to eliminate the boundary effects that 
would otherwise make the decomposition of Q and U into 
E and B ambiguous. 

In Fig. 121 we show the band powers extracted from in- 
dividual runs, using 13 bins logarithmically spaced from 
I — 100 to I = 2000. If the covariance were Gaussian, 
then each estimated band powers would be uncorrelated 
from bin to bin, with rms deviation given by the shaded 
band. The non-Gaussian covariance manifests itself here 
as a positive correlation from bin to bin, and a higher 
variance than one would expect from Gaussian statistics 
alone. This behavior is the Fourier analogue of having 
all of the arcminute scale B modes fluctuate jointly in 
amplitude with the degree scale lenses. 

III. SAMPLE COVARIANCE 
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FIG. 3: Bandpower variance degradation factor Di as a func- 
tion of bandwidth for various choices of the central U of the 
band. The degradation increases with the bandwidth since 
the Gaussian contribution decreases with the number of wave- 
modes. The trend illustrates the effect of the non-Gaussian 
correlation. 



Since the potential fluctuations are Gaussian on the 
large scales that are responsible for lensing as are the 
unlensed i?-modes by assumption, the sample covariance 
can be accurately quantified analytically. 

Consider an ideal, noise-free estimator of the band- 
power 
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dn-B{l)B*{l), 
Aa^ Ji^i 27r 

where A is the survey area in steradians and 



(10) 



(11) 



lei 



is the ?-space area of the band. For a flat spectrum 



A? 



2tt 



const. . 



it is an unbiased estimator of the amplitude 



(12) 



(13) 



where we have used the relationship (27r)^(5(0) = ^ for a 
finite patch of sky. In the limit of bandwidths approach- 
ing = 1, the relationship holds for any underlying 
power spectrum. The band power weight of P/2Tr in 
Eqn. (|10|l is appropriate near the peak at I ~ 10'^. For 
the I <^ 10^ regime where the band power is steeply ris- 
ing but the power is nearly flat (see Fig.^, one can drop 
the weight factors here and still use wide bands. 

Under the assumption of non-overlapping bands in I, 
the sample covariance of the estimator then follows as 



N 



(14) 



The first piece is the Gaussian, or more properly the un- 
connected, contribution 



(5 _ 2(27r)2 j- f If 



(15) 



In the limit of narrow bands and high Z 3> 1, the Gaussian 
piece takes a form that is familiar from Fisher matrix 
studies 



{Ih + l)AZ,/sky 



(A?)^ 



(16) 

where /sky = A /Air is the fraction of sky covered. The 
Gaussian errors mainly reflect a mode counting argu- 
ment. Since the fundamental mode 2tt/A^^^'^ sets the 
spacing of modes in Z-space, Aaj /J27t)'^ is the total num- 
ber of modes in the survey area |l4| . 

The non-Gaussian or connected piece increases the 
band variances and correlates them 



■3 Jhei 



(2^)2 



(27r)2 
where 

bl,^ = W'ih,L)W^i\„L){Cpfctf_^Ct;_^, (17) 
cf;,^ = 14^(1,, l,-L)I^(-l„l,-L)W^(l„l,-L) 



xW^(-l„l,-L)C,t^LQ,-LCrq 
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Note that the both the Gaussian and non-Gaussian pieces 
are the same order in the Cf"^ and Cf"^ power spectra. 

The non-Gaussianity is conveniently quantified by the 
variance degradation factor 



(18) 



which gives the diagonals of the covariance matrix, and 
the correlation matrix 



Ri 



\/ SiiS. 



(19) 



In Tab. 1 we show these quantities as calculated from the 
analytic expression (upper triangle) and 10^ Monte-Carlo 
simulations (lower triangle). With 10^ iterations, the el- 
ements of Rij have converged to a level ranging from 
around 0.01 in the lower bins, to around 0.001 in the 
higher bins (owing to the larger number of modes). The 
remaining discrepancies between the analytic and Monte 
Carlo results are at the 0.01 level and mainly reflect the 
use of the gradient approximation in Eqn. jSJ when de- 
riving the analytic results; the Monte Carlo simulations 
resample Q and U as in Eqn. The close agreement 
between the two provides confirmation that the gradient 
approximation is accurate when computing covariances 
in the range 100 < / < 2000. For these wide bands of 
Al/l ^ 0.25, the non-Gaussian contribution triples the 
variance near the peak of the B power and correlates 
neighboring bands by 70%. 

As the bands widen, the non-Gaussian effects superfi- 
cially appear larger in the covariance matrix. In Fig. O 
we show the degradation factor Di as a function of the 
bandwidth. The nearly linear scaling at high Di can 
be understood as reflecting the linear decrease in the 
Gaussian errors with bandwidth in Eqn. revealing 
a non-Gaussian floor to the variance. For narrow bands, 
this effect is hidden in the band correlations. The net 
effect on the measurements is the same however: when 
combining narrow band measurements to constrain the 
parameters underlying the B power spectrum, their non- 
independence leads to a large degradation in the con- 
straints as we show in the next section. 



IV. SCIENTIFIC IMPACT 

The sample covariance of band powers provides the 
ultimate limitations for measuring the B mode power. 
However until experiments map the B field at a high 
signal-to-noise ratio, detector noise will bring the errors 
closer to Gaussian. 

For a given experiment, the covariance matrix of the 
band powers will include contributions from detector 
noise and the instrumental beam. The sample covariance 
matrix of Eqn. H14|l is replaced with the full covariance 
matrix 



Ci- 



N 



((A^-A^)(A^-A^)} = C-4-5- 



(20) 



Let us make the usual approximation that that noise is 
white and Gaussian. Then llJ 



2(27r)2 



ZTT 



BB 



Ni) 



„;(i + l)e|wHM/81n2 



(21) 



where Ap is the polarization noise variance in a steradian 
and 9 is the FWHM of the beam in radians. 

The impact of the non-Gaussian errors in the power 
spectrum on the errors in a set of cosmological parame- 
ters Pfj^ can be estimated via the Fisher matrix 



(22) 



The errors on a given parameter <j{p^) ~ (F^^)^t{f . Be- 
cause constraints on the cosmological parameters of in- 
terest, e.g. the neutrino mass and the dark energy prop- 
erties, will be affected by constraints from the tempera- 
ture and E power spectra, we defer a full forecast to a 
future work. Instead we calculate the net degradation in 
the signal to noise or equivalently the degradation in the 
errors of the amplitude of the B power spectrum given a 
fixed shape 



C, 



BB 



AC, 



BB 



fid 



(23) 



For reference, the B power spectrum at the peak / ~ 10'^ 
has a sensitivity to the dark energy equation of state w 
of dX/dw\i-i^^ w —0.6 and to the neutrino mass in eV 
of dX/drrii, « —0.5. However these sensitivities should 
not be used for error propagation since they carry an 
/-dependence and depend on what is being fixed by 
the other spectra. For example, the sensitivity to w 
at fixed angular diameter distance to last scattering is 
dX/dw\DA ~ —0.2. Nonetheless the errors on A serve as 
a useful quantification of the overall effect of the non- 
Gaussianity. 

From Eqn. (|22|l the variance in the amplitude becomes 



(24) 



Note that in the Gaussian limit and Al — 1, the variance 
becomes 



a\X) 



/ , 9 /sky 



I - ^-^i 
In Fig. 0] we show the degradation factor 

a\X) 



D. = 



a^{X) 



G 



(25) 



(26) 



as a function of the noise Ap in this Al = \ limit. The 
degradation exceeds a factor of 2 for experiments with 



5 



k 


112 


141 


178 


224 


281 


354 


446 


561 


707 


889 


1120 


1409 


1775 




112 


1.00 


0.03 





04 


0.05 


0.05 


0.06 


0.06 


0.06 


0.07 


0.07 


0.07 


0.08 


0.08 


1 


03 


(1 


04) 


141 


(0.03) 


1.00 





05 


0.06 


0.07 


0.07 


0.08 


0.08 


0.08 


0.09 


0.09 


0.10 


0.10 


1 


04 


(1 


06) 


178 


(0.04) 


(0.05) 


1 


00 


0.09 


0.10 


0.09 


0.10 


0.11 


0.11 


0.12 


0.13 


0.13 


0.13 


1 


07 


(1 


10) 


224 


(0.05) 


(0.07) 


(0 


08) 


1.00 


0.13 


0.13 


0.14 


0.15 


0.15 


0.17 


0.17 


0.18 


0.18 


1 


13 


(1 


15) 


281 


(0.06) 


(0.08) 


(0 


10) 


(0.13) 


1.00 


0.18 


0.18 


0.20 


0.21 


0.23 


0.23 


0.24 


0.24 


1 


21 


(1 


23) 


354 


(0.05) 


(0.08) 


(0 


09) 


(0.13) 


(0.19) 


1.00 


0.22 


0.24 


0.27 


0.27 


0.28 


0.29 


0.29 


1 


28 


(1 


29) 


446 


(0.06) 


(0.08) 


(0 


10) 


(0.14) 


(0.19) 


(0.22) 


1.00 


0.29 


0.30 


0.32 


0.32 


0.33 


0.32 


1 


33 


(1 


35) 


561 


(0.06) 


(0.09) 


(0 


11) 


(0.15) 


(0.21) 


(0.24) 


(0.29) 


1.00 


0.37 


0.40 


0.40 


0.40 


0.39 


1 


55 


(1 


57) 


707 


(0.07) 


(0.09) 


(0 


11) 


(0.15) 


(0.22) 


(0.27) 


(0.30) 


(0.37) 


1.00 


0.47 


0.48 


0.48 


0.47 


1 


76 


(1 


78) 


889 


(0.07) 


(0.09) 


(0 


13) 


(0.17) 


(0.23) 


(0.28) 


(0.32) 


(0.40) 


(0.47) 


1.00 


0.56 


0.56 


0.55 


2 


17 


(2 


17) 


1120 


(0.07) 


(0.10) 


(0 


13) 


(0.17) 


(0.24) 


(0.28) 


(0.32) 


(0.40) 


(0.47) 


(0.55) 


1.00 


0.62 


0.61 


2 


56 


(2 


55) 


1409 


(0.08) 


(0.10) 


(0 


13) 


(0.17) 


(0.24) 


(0.29) 


(0.32) 


(0.40) 


(0.48) 


(0.56) 


(0.62) 


1.00 


0.66 


2 


98 


(2 


94) 


1775 


(0.08) 


(0.10) 


(0 


13) 


(0.17) 


(0.23) 


(0.29) 


(0.32) 


(0.40) 


(0.46) 


(0.54) 


(0.60) (0.66) 


1.00 


3 


19 


(3 


17) 



TABLE I: Correlation matrix Rij and variance degradation factor Di. Parenthetical values are computed from 10^ Monte- 
Carlo simulations and compare well with the analytic calculation. The bands are chosen to be logarithmically spaced and 
non-overlapping with Zmini = 0.89li and Zmaxi = ^mini-i. 




Ap (//K— arcmin) 

FIG. 4: Non- Gaussian degradation on the amplitude A of the 
B power spectrum as a function of detector noise for several 
choices of the FWHM beam Sfwhm. For beams that resolve 
the B power out to the maximum Zmax = 2000 (Sfwhm ^3'), 
the degradation rapidly rises below Ap — 4/iK-arcmin to the 
high signal-to-noise, sample variance limit of ~ 10. Increasing 
'max would further increase the importance of sample vari- 
ance. 



Ap < 4/iK-arcniin and a beam < 3'. At this level, 
the noise within the beam is comparable to the rms of 
the B field [see Eqn. In the limit of zero noise, the 
degradation is a factor of ~ 10 for bands out to I = 2000. 

Note that the degradation factor is insensitive to the 
choice of banding. Computing the variance with the wide 
bands Al/l « 0.25 bands of Tab. 1 causes a negligible 
change in comparison to the accuracy of the underlying 
gradient approximation [see Eqn. |SJ|]. The degradation 
factor is also insensitive to the area or /sky for a large con- 



tiguous region. The degradation can be reduced some- 
what by composing the total area out of many smaller 
patches separated by many degrees on the sky. However 
this strategy will compromise the separability of E and B 
modes. Optimizing a scan strategy against realistic cor- 
related noise, non-Gaussian sample variance and mode 
leakage will require experiment-specific Monte-Carlo sim- 
ulations. 



V. DISCUSSION 

The non-Gaussianity of the B modes in the lensed 
CMB polarization substantially degrades the amount of 
information contained in the B mode power spectrum. 
It both increases the variance of band powers and makes 
them strongly covary across a wide range in I surrounding 
the peak power. Ultimately it will increase the variance 
of the amplitude of the power spectrum by an order of 
magnitude. 

As experiments move from the upper limit and first 
detection stage to using the B mode power spectrum to 
constrain the properties of dark components such as the 
neutrinos and dark energy, this non-Gaussianity will have 
to be included in the analysis. By quantifying the sample 
covariance, we have provided the analytic and numerical 
tools that will be the basis for such an analysis. The 
advantage of the Monte Carlo approach is that it can be 
straightforwardly applied to any estimator of B power. 

In principle, one can include the sample covariance 
in an effective as is done for power spectrum errors 
of the temperature field (e.g. [13 )• However since the 
computation of the covariance is much more costly than 
the computation of the B power spectrum, minimization 
in a large-dimensional cosmological parameter space is 
impractical even with Monte Carlo Markov Chain tech- 
niques. Since most of the parameters affecting the high 
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redshift universe will be fixed prior to these measure- 
ments from the T and E mode spectra, as a first order 
correction one can follow the Fisher matrix approach and 
calculate the effect in a fiducial model. More specifically, 
given the correlation matrix {Rij) and the relative vari- 
ance degradation [Di) in a fiducial model, one can scale 
the covariance matrix to the model B mode spectrum 
(CP^) as calculated from Boltzmann codes. This ap- 
proach would capture the main scaling of the covariance 
through the amplitude of the lensing power spectrum. 
Implementing such a pipeline though is beyond the scope 
of this paper. 

As experiments move from parameter constraints 



based on power spectra to mapping the lensing potential 
001, the non-Gaussianity of the polarization becomes 
the signal and not the noise. The extent to which this 
ultimate goal will be achievable instrumentally and in 
the presence of foregrounds jl6(| awaits the results of the 
upcoming generation of experiments. 
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